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The  invention  of  a direct  method  of  obtaining  candlepower  has  led  to  the  develop- 
ment of  a matiematical  model  of  a simplified  system  for  field  testing  illuminant  items. 

In  addition  to  obtaining  the  light  intensity  independent  of  the  source  distance,  an  in- 
novation for  u sees  sing  the  atmospheric  transmissivity  is  described.  The  light  trans- 
mission is  measured  without  employing  special  equipment  other  than  that  used  for  the 

basic  candlepower  measurements.  The  model  configuration  constraints  are  based  on-  - 
jiregent  field  test  practice.  
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Cont'd 

This  technique  will  allow  two  or  three  people  to  set  up  and  man  a complete  field 
test  system . Test  results  which  formerly  required  extensive  data  reduction  and 
processing  will  be  available  on-site  for  instant  readout.  Besides  increasing  the 
reliability  of  data  and  on-the-spot  readout,  this  system  should  reduce  costs  and 
generally  enhance  illuminant  testing 
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SUMMARY 

The  system  comprises  a complete  method  lor  conducting  field  tests 
on  illumination  flares.  The  physical  model  (fig.  8)  was  selected  for  the 
general  test  requirements  of  artillery  and  aircraft  flare  tests.  It  is  not 
specifically  tailored  to  any  individual  test  requirements.  The  physical 
test  configuration  can  be  optimized  for  specific  test  requirements  while 
coveiing  a specified  flare  drop  area.  In  such  cases,  computer  programs 
can  be  used  to  assess  the  accuracy  of  the  configuration. 

The  cell's  photometer  reproducibility  is  a critical  factor  in  obtaining 
accurate  candlepower  measurements . (Appendix  B has  the  computer  pro- 
gram for  evaluating  transmissivity  error  ve.  sus  cell  reproducibility  for 
the  configuration  shown  in  figure  6 . ) Also  Uie  accuracy  of  obtaining 
candlepower  is  limited  by  the  maximum  fie)  j-of- view . The  angle  subtended 
by  the  cell's  field-of-view  should  be  kept  o a minimum  within  test  con- 
straints . See  appendix  C for  computer  program  to  plot  system  error  versus 
distance  to  the  most  remote  cell . 

The  illumination  falling  on  the  cell  is  measured  to  the  accuracy  of 
the  cell  photometer  itself,  which  will  be  as  good  as  the  system  calibration. 
The  accuracy  of  reading  the  candlepower  will  have  a greater  error  because 
of  the  added  effects  of  changes  with  distance  and  effective  cell  area  vs 
angle.  Candlepower  reading  will  be  limited  to  10  to  20  percent  for  a prac- 
tical system.  Transmissivity  measurements  with  a 2%  photometer  repro- 
ducibility will  have  less  than  10% error. 

There  is  also  error  associated  with  background  illumination.  Error 
caused  by  background  illumination  can  be  largely  accounted  for  by  mea- 
suring and  subtracting  from  the  test,  illumination  levels  because  the  flare 
subtends  only  a small  part  of  the  field-of-view  angle  of  the  cell . Even 
so,  it  is  advisable  to  use  a location  and  configuration  having  small  illum- 
ination levels  compared  to  the  signal . 
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BACKGROUND 

Candlepower  measurements  are  not  made  directly,  but  are 
obtained  by  measuring  the  illumination  and  calculating  the  candlepower 
using  the  inverse  square  law . This  requires  that  the  distance  between 
source  and  sensor  be  known  with  a degree  accuracy.  This  has  been  the 
most  difficult  problem  in  the  field  testing  of  large  flares  because  the  para- 
chute flares  must  be  tracked  as  they  are  deployed  and  descend,  drifting 
with  the  wind . 

The  results  of  tests  of  artillery-  and  aircraft-deployed  illumination 
flares  are  reproducible  when  testing  is  conducted  under  controlled  condi- 
tions. Test  conditions  are  controlled  by  using  tunnels  or  towers  which 
are  available  from  the  tri -services  and  some  contractors.  These  facilities 
are  employed  primarily  m developing  new  systems.  The  typical  informa- 
tion derived  through  testing  is  the  flare's  efficiency,  intensity,  burning 
time  and  spectral  distri  ution  of  light  in  relation  to  the  various  factors 
affecting  the  intrinsic  characteristics  of  the  item  as  a light  source.  When 
deployed  by  aircraft  or  artillery , a descending  flare  produces  spatial  and 
time  variations  in  ground  area  illumination . 

One  facility  for  conducting  field  testing  is  the  Pyrotechnic  Evaluation 
Range  (PER)  located  at  Yuma  Proving  Grounds  (YPG) . The  PER  records 
ground  area  illumination  to  preset  levels,  foot-candles.  From  this  infor- 
mation , a computer  is  programmed  to  produce  all  required  test  data . 

A major  problem  encountered  in  field  testing  is  controlling  or  ac- 
counting for  environmental  factors  affecting  the  test.  Previous  methods 
employed  to  overcome  this  problem  have  required  many  people  to  trans- 
port. locate,  and  man  the  equipment.  The  quality  of  the  test  data,  which 
are  often  incomplete  or  unreliable,  does  not  justify  the  cost  of  this  method. 
For  example,  evaluating  the  flare's  candlepower  in  the  field  requires  that 
the  flare's  position  be  determined  versus  time.  The  PER  uses  three  CINE 
theodolites  for  this  purpose,  and  many  people  to  operate  them.  Another 
field  testing  problem  is  correlating  times  and  coordinating  data  systems 
with  the  firing  of  the  flare.  During  field  tests  performed  at  Ft.  Hood, 
manned  photometers  were  deployed  over  large  areas;  stations  were  separ- 
ated by  thousands  of  meters. 


New,  ultrasensitive,  prototype  photometers  were  used  with  L-Band 
RF  Links  to  record  remote  location  information  in  a telemetry  van . The 
flares  were  tracked  by  optical  techniques.  Typically,  some  data  is  lost 
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during  tests  of  this  kind  because  the  position  of  the  item  is  lost  temporarily 
during  its  burning.  Also,  fibres  descending  and  traverse  with  the  wind 
occasionally  saturate  near  st»  ions  and  produce  weak  signals  at  distant 
stations. 

To  evaluate  the  light  source  intensity  with  any  degree  of  accuracy, 
the  atmospheric  transmissivity  of  the  test  range  environment  must  be 
assessed.  In  the  past,  weather  reports  provided  crude  estimates  of  visi- 
bility limitations.  Also,  since  the  varying  configuration  of  a flame  pro- 
duces a flickering,  nonisotropic  light,  source  intensity  can  be  determined 
reliably  only  when  it  is  derived  through  spatial  integration  of  candlepower, 
i.e. , average  light  intensity  when  viewed  from  all  directions. 

The  method  described  in  this  report  uses  an  array  of  highly  sensi- 
tive photometers  incorporating  automatic  change  o7  scale  circuitry  and 
automatic  calibration  to  operate  without  a man  in  the  loop. 

DEVELOPMENT 

System  development  involves  two  innovations  to  the  state-of-the-art 
of  flare  testing.  The  first  is  a method  of  determining  candlepower  of  tne 
flare  item  without  knowing  its  position  or  distance . The  second  is  a means 
of  assessing  the  transmissivity  of  the  atmosphere  over  the  test  area  to  an 
accuracy  better  than  10%.  The  candlepower  measurement  is  corrected  by 
the  transmission  factor  to  yie'd  the  true  field  candlepower  of  the  flare.  The 
photometers  record  ground  illumination . All  other  parameters  are  derived 
through  data  processing.  These  measurements  are  taken  in  real  time  as 
the  item  burns  and  descends  by  parachu'3. 

The  photometers  have  a fixed  field-of-view  and  a position  that  is 
determined  by  the  system  equations  and  mathematical  model.  A restricted 
field-of  -view  and  optimum  angle-of-iook  are  required  to  reduce  system  errors 
caused  by  background  illumination  and  field  test  configuration.  The  field- 
of  view  for  each  cell  will  be  designed  through  the  use  of  a series  of  baffles 
to  view  only  the  region  of  operation  of  the  item.  The  effects  of  background 
illumination  will  be  further  reduced  by  an  automatic  calibration  and  cor- 
rection system. 

Calibrating  of  the  system  is  a three-step  procedure.  First,  the 
optical  equipment  is  calibrated  with  a laboratory  standard.  Second,  equip- 
ment in  the  field  is  calibrated  with  an  optical  transfer  standard.  Third,  the 
signal-conditioning  electronics  and  RF  link  are  electronically  calibrated . 


The  field  photometers  consist  of  an  ICI-Y  photocell  which  is  corrected 
to  ofc.jin  spectral  response  equal  to  the  average  human  eye,  a filter,  a 
cone,  and  electronics  (fig.  11 . A silicon  photovoltaic  barrier  layer  photo- 
cell is  used  for  its  stability  and  small  temperature  coefficient.  The  cell 
drives  a low-noise  operational  amplifier  which  is  used  in  its  inverting 
mode  to  enhance  the  system's  dynamic  charge.  Operating  the  coll  into  an 
effective  zero  impedance  extends  the  cell's  linear  range.  Three  operational 
amplifiers  (G1  ,%ru  G3,  figure  1)  measure  resolution  and  dynamic  range. 
The  automatic  circuitry  selects  the  appropriate  amplifier  and  outputs  data 
to  a co--  ,-r Jonal  telemetry  system.  The  photometer  output  is  transmitted 
via  an  L i ind  PCM/FM/FM  RF  link  from  each  rer  ite  field  sensor  to  a cen- 
tal recording  van  where  the  data  are  processed  and  recorded . Tc  improve 
system  accuracy  each  photometer  has  a pulse  code  modulator  (PCM)  which 
digitizes  be^re  transmission.  The  PCM  word  length  is  set  to  reduce  quan- 
tization errors  to  meet  system  requirements.  Parity  is  used  to  detect  trans- 
mission link  errors.  E_ch  photometer  incorporates  a voltage  controlled 
oscillator  (VCO)  which  contains  photometer  scale  information.  An  auto- 
m-  ‘>c  change  of  seal?  is  used  to  improve  the  system  dynamic  range.  The 
outi?u*s  of  the  PCM  anc  VCO  modulate  the  L-band  RF  transmitter. 

. automatic  circuitry  selects  the  appropriate  system  gain  to  optim- 
ize •’ ..  jtput  data,  provides  scale  information  on  a separate  VCO  channel, 
ar  ourcs  data  calibration  by  an  automatic  calibration  interrupt. 

The  transmitted  signals  from  each  field  photometer  are  received  by 
an  L-band  antenna  and  receiver  at  the  central  receiving  station  Ifig . 2) . 
The  receiver  output  is  applied  to  two  discriminators . One  discriminator 
cor.  = jins  the  illumination  data,  in  PCM  format,  and  the  other  provides  the 
automatic  range  information. 

The  processor  operates  on  both  the  data  and  scale  information  from 
each  remote  photometer  and  outputs  to  a visual  display  the  real  time  illum- 
ination levels  (E)  for  each  remote  site. 

A simplified  general  field  test  configuration , using  three  remote 
photometers  and  a source  flare,  is  illustrated  in  figure  3.  Dj , D2 , and  D, 
represent  the  distance  from  the  source  flare  to  each  respective  photometer 
with  E1 . Ej  , and  Ej  representing  the  levels  of  illuuinadnn  they  receive . 
These  illumination  levels  are  transmitted  to  the  central  recording  site  and 
impressed  on  the  input  of  the  processor . At  the  beginning  and  at  the  end 
of  the  test,  the  processor  calculates  tne  transmissivity  (T)  from  E, , E2 , 
and  Ej . Using  the  transmissivity  and  any  two  of  the  illumination  levels, 
the  processor  calculates  the  intensity  of  the  source  (I) . 


Candlepower  Measurements 

Candlepower  is  calculated  from  the  difference  of  two  photometer 
readings.  A typical  geometry  is  shown  in  figure  4.  Photometers  are  lo- 
cated at  fixed  positions  E2  and  Es  with  a fixed  field-of-view  and  directed 
as  shown. 

This  system  uses  an  in-line  array  of  two  photometers  to  make  the 
flare  'tern  appear  as  an  isotropic  radiator.  Two  photometers  are  placed  in 
a line  at  a distance  from  the  anticipated  position  of  the  flare's  trajectory. 
Given  these  conditions,  the  candlepower  of  the  test  item  can  be  determined 
by  a relation  of  the  two  illumination  levels  recorded  by  the  two  photometers 
E2  and  Ej  and  the  fixed  distance  between  them.  This  eliminates  the  need 
to  track  the  flare  and  measures  its  distance , with  respect  to  time . from 
the  photometers. 

Consider  the  simple  case  (tig.  5}  where  the  fare  (source)  is  at 
position  0.  The  expression  for  the  illumination  being  seen  by  the  photo- 
meters is: 


E.  = -. 


E*  =1 


D,  = D,  + d 


To  solve  fo;  the  candlepower  (I)  of  the  '.'are,  the  following  equations 
are  used  (distances  Itj  and  D2  are  unknown)- 


from  equation  i 


from  equation  2 


where 


I = E,D* 


I = E,Df 


E = Iliuminaticr 


I - intensity  of  source  in  candlepower 
D = Distance. 


jwto?'??;  - - - ►*■- 
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Then 

E,D|  - E2D|  = 0. 

Using  equation  3 to  eliminate  Ds . 

Es  (Dj  + d)*  - Ej D*  = 0. 
Solve  for  D2  using  quadratic 


-d  (Ej  +J  E2Es ) 
°1=  (E,  -E,) 


(6) 


(7) 


(8) 


Substitute  D2  to  solve  for  1 


I =■  E3  (D,  + d)1 


I = F, 


d d_(Ej  4 Ve2  E,) 

e,-e2 


(9) 


(10) 


Put  in  simplest  form 


I = 


(11) 


Therefore,  knowing  the  illumination  levels  E2  and  Es  and  the  fixed  distance, 
d, separating  photometers  #2  and  #3,  it  is  possible  to  compute  candlepower  I. 


This  simplified  case  does  not  represent  a typical  flare  field  test 
configuration  because  the  flare  in  the  field  test  is  not  directly  in  line  with 
the  photometers-  Flares  are  deployed  to  burn  from  a height  of  approxi- 
mately 2,000  feet  to  800  feet  and  then  extinguish.  Measurement  accuracy  is 
improved  by  aiming  the  cells  to  view  this  region  of  operation  (fig.  4)  op- 
timizing the  angle  at  which  the  light  rs/s  impinge  on  the  cell.  This  angle 
is  critical  because  if  it  is  not  optimum , it  will  reduce  the  effective  area  of 
the  cell  which  in  turn  reduces  the  illumination  level  that  the  ceil  experiences. 
A light  incidence  angle  of  10°  off  ceil  axis  normal  would  cause  an  error  as 
high  as  4%.  To  partially  compensate  for  this  error,  the  cells  will  be  aimed 
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at  a point  between  the  800  ft  and  2,000  ft  elevation.  For  the  test  configur- 
ation of  figure  4,  cells  Ej  and  Es  are  set  9,93°  and  7.25°  respectively, 
above  horizontal. 

The  computer  program  in  appendix  A computes  the  expected  accuracy 
of  the  system  and  the  optical  separation  between  stations  E2  and  E3  for  the 
configuration  inputed.  For  the  test  configuration  of  figu;e  4.  the  optical 
separation,  d,  is  computed  using  equation  (3): 

d - D3  - D2  (12) 

where  Dj  - 11088.73  ft 
and  Dj=  8121.58  ft 
then  d = 2967.15  ft 

This  value  of  d applies  only  to  the  configuration  described  in  this  report. 

A new  value  of  d would  be  computed  when  the  configuration  is  changed . 

Besides  eliminating  the  tracking  equipment  and  the  additional  people 
needed  to  set  up  and  man  the  tracking  equipment,  this  method  simplifies 
the  data  recording  and  evaluation.  With  the  old  method,  extensive  data 
processing  was  required;  with  this  system  the  data  are  processed  while 
the  test  is  in  progress.  The  ground  illumination.  E,  and  E3  . is  being  re- 
corded and  plotted  in  real  time.  The  flare's  intensity  (I)  is  also  being 
computed  and  plotted  using  the  equation  (11)  and  E2  . E3  and  the  system 
constant,  d.  This  requires  a data  processor  consisting  of  three  integrated 
circuit  chips  for  generating  the  flare's  intensity  in  real  time  analog  form. 

In  actual  field  tests  more  than  one  array  of  two  photometers  may  be  deployed 
depending  on  test  requirements.  There  are  no  stringent  requirements  on 
the  photometer/ receiving  station  RF  link;  L or  S band  telemetry  trans- 
mitters have  been  used  with  FM/FM  modulation.  Depending  on  system 
accuracy  required,  either  FM/FM  or  PCM/FM  can  be  used. 

Transmissivity  Measurements 

Transmissivity  is  calculated  by  finding  the  difference  between  iJ' ani- 
mation measurements.  The  difference  between  the  two  measurements  in 
the  array  from  the  same  source  indicates  the  transmission  of  light  on  the 
optical  path. 
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The  system  for  measuring  transmissivity  uses  highly  sensitive 
photometers  with  automatic  change-of-scalo  circuitry  to  extend  the  dynamic 
range  without  limiting  instrumental  accuracy.  The  in-line  array  of  three 
photometers  (fig.  6)  yields  a method  of  assessing  the  loss  of  light  caused 
by  atmospheric  attenuation  to  within  10%.  If  a large  flare  is  burned  as 
shown  in  figure  6 before  and  after  each  test,  the  transmissivity  of  the  test 
range  environment  can  be  measured.  Figure  7 shows  a typical  test  config- 
uration . 


The  general  expression  for  determining  the  intensity  (I)  of  a flare 
with  an  in-line  array  of  two  photometers  is 


A 


r = d/Cl^lEi  - l/^EJ 


(13) 


where  d is  the  known  optical  separation  between  sensors  II  and  I , and  E2 
and  Ej  are  the  illumination  levels  (fig.  6) . This  general  expression  can 
be  modified  to  include  atmospheric  attenuation  by  scaling  the  illumination 
levels  as  a function  of  T (transmissivity)  and  distance  from  source. 


I5  = d/^1/  V(E,/TD  - 1 \l  Ej/T0"^ 


(14) 


A similar  equation  can  be  obtained  using  the  illumination  data  from 
sensors  II  and  III . 


I * — d ( 1 / 1|(E,/Td+D  -1/  \/  E2/TD 


\[e,/TD  ) 


(15) 


Eliminating  I in  equations  (14)  and  (15)  yields 


(i/||"e:/TD  -1/  if  E,/TD‘d  ) 

- (]/f  *>' 


TD+d  -1/  « Ej/T 


'j  E, 


(16) 


lying 


By  multiplying  numerator  and  denominator  by  1/  ]j  TD  and  simpli- 


,d/2 


.-d/2 


P7  ^ 


(17) 


8 


Solve  for  T by  changing  variables  and  using  quadratic  equation. 


-- U|«-j/e,  ♦ Ej  -\L 


The  solution  of  T yields  the  transmissivity  for  the  test  range  en- 
vironment, which  can  be  inserted  in  equation  (14)  or  (15)  along  with  the 
measured  values  of  E, , Ej , and  the  fixed  value  of  d to  obtain  a more  ac- 
curate field  measurement  of  the  flare  candlepower.  This  system  assesses 
the  atmospheric  transmissivity  of  the  test  range  before  a test  is  conducted 
Test  data  are  accurate  to  within  10%.  Transmissivity  tests  are  performed 
with  the  equipment  provided  for  the  basic  illumination  measurements,  so 
the  system  is  highly  efficient  and  very  economical. 
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Figure  2.  Receiving  station  block  diagram  (single  channel) 


configuration  for  flare  measurement 


e 6.  General  test  configuration 


Typical  test  configuration 


•^i  K-^ 


1 r.E-i  — appendix  a 

2 nCI *1  TL£  - PI-OCftA.l  TO  FlliL  "CPST  CASE  f:  * 07  IN 

3 nE;.j CAfJELEPO*.  EP  CALCLt.ATJO.i^.IEuLECTII.C 

„ r-ai TRANSMISSIVITY  FI  D.J  TEST  configuration 

<■  TFJi AS  SKO'.-l  1.1  Flo.  ‘4- 

10  PRINT  "***•***•«****•*•*•••***••«»••*•»*»***••*••*—•«•**«* 

20  PRINT 
30  PRINT 
40  PRINT 

45  PRINT  "PRESS  RETURN  TO  CONTINUE  PROGRAM" 

46  INPUT  A 

47  IF  A*»0  THEN  570 

50  PRINT  "INTENSITY  OF  FLARE”? 

51  INPUT  I 

52  PRINT  "MIDDLE  OF  CUBF.  CX-C00RDIMATE1 " 

53  INPUT  XI 

60  PRINT  "ERROR  OF  CELLS  IN  DECIMAL  FORM"! 

61  INPUT  E 

62  PRINT  "CELL  SEPARATION"? 

63  INPUT  D 

65  PRINT  "PRESS  0 TO  CHANGE  VARIASLLS-  OTHER VI SE  PRESS  RETURN" 

67  INPUT  J 

66  IF  J<»0  THEN  10 

70  PRINT  "DISTANCE  IN  X DIRECTION  FROM  FURTHEST  CELLS"? 

71  input 

72  PRINT  "Y»"? 

73  INPUT  Y 

74  PRINT  "Z*”? 

75  INPUT  Z 

80  LET  DS-SQR<X1-XI»!  -96000£*0  6>-SCR(  (XI  - D>  »<X!  -D>  » I -960D0E+0  6) 

81  LET  T1*ATN«1400/(X1-D>) 

82  LET  T2-ATN<1400/X1> 

130  GOSUB  480 

215  LET  B-D 

216  LET  H3«H! 

260  PRINT  "INTENDED  SEPARATION  ’D'  *"D5 
270  'PRINT  "INTENSITY  OF  FLARE-1"! 

280  PRINT  “CELL  ERROP.*"ABS(  £*  10E > 

300  PRINT  "DISTANCE  BETWEEN  C£LLS-”B 

310  PRINT  "ERROR  OF  SYSTEH*"H3 

320  PRINT  "LIGHT  INTENSITY  SYSTEM  FIND$-"P 

360  PRINT  " — - “ 

370  GOTO  65 


400  REM D- SEPARATION 

402  REM D5-1NTENDED  SEPARATION 

REM X-Y-Z-COORDI  NATES 

410  REM 1 ■ INTERS!  TY 

415  REM E- ERROR  OF  CELLS 

420  REM SUBROUTINE  TO  CALCULATE  ERROR 


421  LET  VI-SaR<X-X-Z-Z> 

422  LET  V2*ATN(Z/X) 

423  LET  V3-V2-T2 

424  LET  V4-VI-C0S(V3> 

425  REM HI -ERROR  OF  SYSTEM 

426  LET  V6-SQR< (X-D5-<X-CJ*Z-Z> 

427  LET  V7-ATN(Z/  <X-D) ) 

428  LET  V8-V7-T1 

429  LET  V9-V6«C0S< V8> 

430  GOSUB  580 
440  LET  M2-H1 
445  LET  P5-P 
450  LET  E— E 

A4«  J T-'  OC-O 

460  GOSUB  500 

478  IF  HI  >•  K2  THEN  560 

480  LET  H1-H2 

485  LET  P-PS 

490  GOTO  560 

500  LET  DI-<X-D>»(X-D>*Y-Y*Z«Z 

S10  LET  D2-X»X»Y»Y»Z*Z 

520  LET  El-<I/Dl>*<V9/SQR(Dl>>*<t»E> 

530  LET  E2-<I/D2)*(V4/S0R(D2>?»J1-E) 

540  LET  PI*D5/<  l/SOR<  E2)-1/S0R(  El  > > 

545  LET  P-P14P1  ' 

550  let  hi- u?s<(«p-n/n«iee> 

560  RETURN 
570  END 
READY 


( 


PRECEDING  PADS  BLaMK-SDT  FIUCD 


00M 


21 


TT  r*-; 


y ; =?; 


RUN 

*********************** ************************* ******** 


PRESS  RETURN  TO  CONTINUE  PROGRAM 
7 

INTENSI TY  OF  FLARE?  J 000000 
MIDDLE  OF  CUBE  (X- COORDINATE) 

?1 1000 

ERROR  OF  CELLS  IN  DECIMAL  FORM?. 02 
CELL  SEPARATION? 3000 

PRESS  9 TO  CHANGE  VARI ABL ES* 0 TH ERVI S E PRESS-  RETURN 
7 

DISTANCE  :N  X DIRECTION  FROM  FURTHEST  CELLS?! 1000 

Y=?0 

Z=? 1 400 

INTENDED  SEPARATION  'D'  = 2967.  157 
INTENSITY  OF  FLARE=  1000000 
CELL  ERROR-  ?. 

DISTANCE  BETWEEN  CELLS=  3000 

ERROR  OF  SYSTE4=  14.29149 

LIGHT  iNTENSI  TY  SYSTEM  FINDS*  1142915 


PRESS  9 TO  CHANGE  VAP.I  ABL  E S>  0 TH  LRU  I S E PRESS  RETURN 
?9 

******************************************************** 


PRESS  RETURN  TO  CONTINUE  PROGRAM 
?9 

STOP  AT  LINE  570 
READY 


wyM  < tv*#*?*  rwv^y^,u!iw ggwgpwf* 


1 REM TITLE  - SPECIAL  PROGRAM  FOR  FINDING  WORST  CASE 

2 REM TRANSMISSIVITY  ERROR  VS.  CELL  REPP.O  DUCABLI TY 

3 REM  USING  CONFIGURATION  SHOWN  IN  FIG. 6 

12  PRINT 

PRINT 

PRINT  " TRANSMISSIVITY  PROGRAM  I" 

PRINT 
PRINT 

PRINT  “INPUT  THE  FOLLOWING;  T,H,  E"; 

INPUT  T,H,  E 
IF  T<e  THEM  17 
IF  H<0  THEN  17 
17 


13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

70 

71 

72 
75 
60 
85 
90 
95 
105 
110 


IF  L<0  THEN 
PRINT 


1 — D3—  ! — D2—  l- 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 

print  ,,t»w;t;m#ch.l  ERROR- **JE 

PRINT 
PRINT 
PRINT 
PRINT 
PRINT 

PRINT  “INPUT  C1#>D2,D3" 

INPUT  D! • D2> D3 

111  IF  D2< > D3  THEN  105 

112  IF  D1  < 0 THEN  10S 

113  IF  D3<0  THEN  .05 

114  IF  D1  < >0  THEN  119 
U5IF  D3< >0  THEN  119 
116  GOTO  10 

119  LET  R*0 

120  LET  L3* ( DI  ♦D2^D3) » 2+H-H 
130  LET  L2a ( D2+D1 > *2+H*K 
140  LET  LI-D1*DUH*H 


TEST  CONFIGURATION" 


-Dl- 


; H" FT" 


D3  AND  D2  MUST  BE  EQUAL” 


’’TYPE  0 FOR  D1  j D2*  D3  TO  RESTART  PROGRAM" 


150 

FOR 

X— I TO  1 STEP  C 

160 

FOR 

Y*-l  TO  1 STEP  2 

170 

FOR 

2*-l  TO  1 STEP  2 

160 

LET 

A*E*2 

190 

LET 

8-  £*Y 

20C 

LET 

C*  £*X 

210 

LET 

El  * 1 00C00. /L 1 *T* ( SGRCL 1 >/ 5260  > * < 1 ♦A) 

220 

LET 

E23 1 00000  «/L2*Tt ($CR(L2)/ 528  C > * < 1 ♦E) 

230 

LET 

E3- 1 000  00./ L 3*  T* < SCR<L3)/S260  >*<  I ♦ C > 

240 

LET 

T1  -SQR<  E3/E2) 

250 

LET 

T2-SQRC  E3/E2-SQR<  E3/EI  >) 

260  LET  T3-CTUT2)  «<  2/ <03/5280  >) 
270  LET  R1»(T3-T)/T^100 

IF  ABS<  R1 ) >ABS<  R)  THEN  00 
NEXT  Z 


260 

290 

300 

310 

320 

330 

340 

350 

2$e 

370 

390 

400 

410 

420 


NEXT  Y 
NEXT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
GOTO  9 5 
LET  R-P.l 
GOTO  290 
END 


Di  ■••;  Di ; "D2«-;  D2;  -D3--;  D3 

error  of  t--;r 


2k 


I 

f 


APPENDIX  C 

Program  to  Plot  System  Error  vs  Distance  to  Farthest  Cell 


r 


'*'*’■*  y ,1^.*^ ■*? iw .-AT' x .•V.'^v-? 


1 REM  APPEND iy  C 

2 REM  TITLE  PROGRAM  TO  PLOT  SYSTEM  ERROR  VS 

DISTANCE  TO  FARTHEST  CELL 

LIST 

3 LET  0-1250 
10  LET  J-2000 

IS  LET  E« 100000. 

17  LET  N—2.00000E-02 
20  LET  L-300B 

25  DEF  rNXIX>-<X-L)»2*G»2*J»2 
30  DEE  FNY(X><“Xt2*G*2*  J*2 
35  DEE  ENZ<X)-SOftt<X-L>»2*Jt2) 

40  DEE  mV(X)-ATN<J/<X-L))-ATN<1400/CX-L>> 

45  DEE  FNT<X»-FNECX>»COSCFNWlX>> 

50  DEE  FNHtX>-SjR<X*8M.960e0E<06>-SOR«<X-L>ta*!.96000E*06) 

55  DEE  FNS<X>-30R<X»2*Jt2> 

60  DEE  ENRtX>»ATNt J/XI-ATNC 1400/X) 

65  DEE  FNO<X>-FNS<X>*COS<FNR(r>> 

70  DEF  FNPCX)«<E/FNX<X>>*(FNT(X>/S0R<FNX<X>)>*<1*H> 

75  DEF  FNOCXI-C  E/FNYCX) >*C  FNO(X) /SCRC FNYCXJ  > )•€ l-M> 

00  DEF  FNNIX)«<  FNNCX) /< l/SQRC FN0CX1 >- t /SQRC FMP(X> ) ) ) »8 
100  DEE  FNFCXI-<FNN<X>-E>/E*100 
IS*  DIN  EC  10] 

400  LET  Rl-0 
900  LET  Ll-0 
600  LET  01-0 

700  PRINT  “PLEASE  INPUT  THE  FOLLOWING  PARAMETERS  I " 

000  PRINT  “LEFT  X-ENDPOINT"! 

900  INPUT  A 

1000  PRINT  “RIGHT  X-ENDPOINT"! 

1100  INPUT  B 

1000  PRINT  “X-SPACING"! 

1300  INPUT  D 

1400  PRINT  "THE  NUMBER  OF  UNDEFINED  POINTS  CIF  NONE, ENTER  0>"! 

1500  INPUT  N9 

1600  IE  Nv-0  THEN  8100 

1700  PRINT  “ENTER  THE  UNDEFINED  POINTS,  FOLLOWING  EACH  WITH  A RETURN" 

1 800  FOR  K7-1  TO  N9 
1900  INPUT  ZCK7] 

8000  NEXT  K7 

8100  DEE  FNG<X>«INTt  C Y7-L1 >/Dl 5>* 1 5 
8800  LET  L2-R2-FNMA) 

8300  FOR  X-A  TO  B STEP  D 
8400  FOR  I- 1 TO  N9 
8500  IF  X-ZCI]  THEN  3100 
8600  NEXT  I 

2700  IF  FNF(X)>L8  THEN  2900 
8800  LET  LB-FWFCX) 

2900  IF  FNF(X)«R8  THEN  31J0 
3000  LET  RB-FNF<X> 

3100  NEXT  X 

3200  IF  L8<0  THEN  3500 

3300  LET  R1-R2 

3400  GOTO  3990 

3500  IF  R8>0  THEN  3700 

3600  GOTO  3800 

3700  LET  R1-R8 

3800  LET  L1-L8 

3900  LET  DI-<R1-L|)SS0 

4000  IF  L!«R1  THIN  4300 

4100  PRINT  "THIS  IS  THE  FUNCTION  T-CONSTANT." 

4800  STOP 
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— A I-'*-  -*-  ',--1 


fe  '-t  . .. --R, 


«Mt 

4*W 

MM 

«TM 

«M 

•M 

MM 

SIM 

MM 

MM 

MM 

MM 

MM 

STM 


Miirr 

mimt 

MIMT 

MIMT 


-•Cl  V*-UK  or  TME  ftMCTIOM  1S«JL* 

-21  UiiTIf!  ^UI  0r  ™*  WMCTIOM  IS«Jft» 

JJ**  SPACiNa  o«  the  t-axis  is**idi 


lit  r<iim-Li/Dii.5)tis 
IK  <•  I THEN  S9M 
IE  A/D»8  THEM  !)H 
LET  «l>l 

IE  Ll«*  THEM  33M 
PRINT  TAAC  EJ 1 
EEIHT 

00  TO  78  99 

E0»  1*1  TO  tNT<ASD-.S> 
EEIHT  TA*(E>|-4- 

HOCT  I 


*M 
Ml 
EM* 
• IN 
MM 
UM 


LET  Bl>* 

EOR  X.A  TO  S STEP  0 

IE  rx  ! * MIME-  II  THE*  63M 

IE  AlS<X>»|.MM*E-tS  THEM  till 

LET  X-C 

EEIHT  X. 


MM  EOE  p-l  to  HE 
•»•  IE  X*E'P1  THEM  7SM 
**M  IE  XEE  THEM  7399 
«7M  EOE  1 t-  1 TO  S( 

MM  PEI MT  “4-i 
MEXT  It 
7ME  LET  0-1 
TIM  EEIHT  -T- 
TSM  00TO  MET 
73**  PEIHT  TA*(EII*4- 
741*  0070  9117 

T5M  HIXT  P 

T»M  IE  X«<X4S>>(  THEN  9699 
TIM  IE  X«-0/t  THEM  EM* 

7*M  EOE  t-E  TO  St 
TEM  IE  Bit*  THEN  6EM 
M*»  LET  T7-mr<x> 

•IM*  IE  EHEtXl-l 4| s THEM  SSM 
**••  IE  I*1S*E  THEM  8?t* 

*3**  FEINT  *4-1 
cm  GOTO  8896 
*3**  PEIHT  “4”l 
MM  0070  88M 
*7 EE  PEIHT  -0-1 
8RM  HIXT  I 

•EM  IE  I.ISIE  THEM  EIE* 

•ME  PEIHT  -4-1 

•I**  PEIHT  -7“ 

EtM  LET  0-1 

E3M  IE  (Bfll.l  THEN  E987 

•m  IE  (0141). t THEM  35*3 

*3M  IE  (Sl4|).]  THEN  EEI» 

****  IE  X4(X-O>40  THEM  E8M 

EE**  IE  X «.  D/»  THEM  78V* 

E8EI  LET  T7-EHEC  X> 

•EM  IE  EMO(X) 4E  THEM  9986 
•9*1  IE  EMC  ->.E  THEM  998* 

EE8i  PEIHT  -«*8(  EHG(  X7>  !“•“)  T8B(  E)  I "4— 

•9*3  GOTO  EEE7 

»M*  PEIHT  TUIEII-4- 

•913  GOTO  0987 

•EM  PEIHT  TEE(  El  |— 4— ) TABf  EM0(X  > >1—4— 

•*•7  NEXT  X 

•M*  IE  X 4.  8 THEM  9017 
•*99  |p  -X/D»6  .HEM  9917 
MI*  EOE  1. 1 TO  lNT(-X/0-.S> 

••II  EEIHT  TA8IEU-4- 

Ml*  HIXT  : 

MIS  LET  Sl*t 
MI*  PRINT 
MIS  GOTO  T8M 
Mil  PEIHT  TMI7II-4- 
MIT  PEIHT  TM(EII*X“ 
t*I»  IP  9-1  THIN  9*18 
MI9  ETOP 
MM  PEIHT 
Mil  PEIHT 
MM  PEIHT 
MIS  EOE  TO  39 
•*««  PEIHT  *4-1 


9*9*  PEIHT  T" 

MET  PEIHT 
MM  PEIHT 

mS  mV  " *'"CC  ™I  T*«'S  15  0"  the  GRAPH. 

Mil  END 


MAST 
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